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Abstract 

Mathematical models have recently been used to cast doubt on the biotic origin of 
stromatolites. Here by contrast we propose a biotic model for stromatolite morpho- 
genesis which considers the relationship between upward growth of a phototropic 
or phototactic biofilm (v) and mineral accretion normal to the surface (A). These 
processes are sufficient to account for the growth and form of many ancient stro- 
matolities. Domical stromatolites form when v is less than or comparable to A. 
Coniform structures with thickened apical zones, typical of Conophyton, form when 
V >> X. More angular coniform structures, similar to the stromatolites claimed as 
the oldest macroscopic evidence of life, form when v >>> A. 



1 Introduction 



Stromatolites preserve the only macroscopic evidence of life prior to the ap- 
pearance of macro-algae [1]. The biogenicity of stromatolites older than 3.2 
Ga^ is unclear [2,3,4,5]. If they are indeed biotic, they are the oldest mor- 
phological evidence for life, now that the identification of 3.3 to 3.5 Ga micro- 
fossils [6] has been challenged [7,8]. Here we propose a mathematical model 
for stromatolite morphogenesis that endorses a biotic origin for coniform stro- 
matolites. It analyses interaction between upward growth of a phototropic 

^ 1 Ga = 10^ years. 



Preprint submitted to Elsevier Science 



8 February 2008 



or phototactic microbial mat and mineral accretion normal to the snrface of 
the mat. Domical structures are formed when mineral accretion dominates. 
When vertical growth dominates, coniform structures evolve that reproduce 
the features of Conophyton, a stromatolite that flourished in certain low- 
sedimentation environments for much of the Proterozoic [9]. Increasing the 
dominance of vertical growth produces sharply-peaked conical forms, com- 
parable to coniform stromatolites described from the 3.45 Ga Warrawoona 
Group, Western Australia [5]. 

Some authors prefer to avoid a genetic definition for stromatolites [10], but we 
regard them as laminated microbialites [11], biomechanically [12] and func- 
tionally [13] analogous to lithified sessile organisms, such as colonial corals, in 
which hving tissue is restricted to the surfacel [14]. In stromatohtes the hving 
tissue is a benthic microbial community (BMC) [15]. BMCs range from com- 
munities composed of a single species to complex trophic networks of photoau- 
totrophs, chemoatotrophs and heterotrophs [15] in which species composition 
and diversity may change in response to environmental conditions [16] In cases 
where the BMC includes photoautotrophs, such as cyanobacteria, the stroma- 
tolite form represents a record of that community's response to light. There 
have been very few previous attempts to model stromatolite morphogenesis 
mathematically. Verrecchia [17] proposed a simulation model for microstro- 
matolites in calcrete crusts using a diffusion-limited aggregation model [18], 
but this is only relevant to modelling complexly branching stromatolites. In- 
sights into morphogenesis of simpler forms may be gained using the interface 
evolution equation of Kardar, Parisi and Zhang (KPZ) [19] which contains pa- 
rameters for surface-normal accretion, surface tension, and noise. Grotzinger 
and Rothman [20] attempted to simulate stromatolite form using a modified 
KPZ equation with explicit vertical growth. In their model vertical growth was 
considered to be due to the deposition of suspended sediment, surface-normal 
accretion was due to chemical precipitation, surface tension effects were re- 
lated to both diffusive smoothing of the settled sediment and chemical precip- 
itation, and uncorrelated random noise represented surface heterogeneity and 
environmental fluctuations. Their model simulated the structure of a supposed 
stromatohte from the Cowles Lake Formation, Canada, and this led them to 
conclude that some and perhaps many stromatolites may be accounted for 
exclusively by abiotic processes. However, that model was subsequently mod- 
ified to include a biotic process, mat growth, along with mineral precipitation 
in the surface-normal growth parameter [10]. As both mineral accretion and 
biological growth were linked in their surface-normal growth parameter, this 
model was unable to discriminate biotic effects. A rather different application 
of the KPZ equation has been proposed [21] in which the effects of microbial 
growth are included in the vertical growth parameter of the equation. This has 
been used to simulate the morphogenesis of stromatolites from Marion Lake, 
South Austraha. 
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2 The model 



The biotic interpretation of fossil stromatolites is widely accepted, despite the 
fact that they rarely preserve any remains of the BMC which formed them 
[10]. As a result, attention has been focussed on how biotic stromatolites might 
be distinguished from abiotic accretions such as tufa, speleothems, and cal- 
crete [10]. Several Proterozoic stromatolite forms grew in environments of low- 
sedimentation and their formation seems to have been due to the growth of a 
BMC, containing photosynthetic bacteria, and accretion of calcium carbonate 
in the resulting biofilm [1,3,9]. For forms which lack evidence of detrital mate- 
rial being trapped or bound by the BMC we propose a model for stromatolite 
morphogenesis which involves two processes only: 

(1) upward growth of a phototropic or phototactic BMC, 

(2) mineral accretion normal to the surface. 

The function h{x, t) represents the height of the profile above a horizontal 
baseline which evolves in time t according to the equation 



The co-ordinate x measures the distance along the baseline. It is also equiva- 
lently a radial co-ordinate in the baseplane for circularly symmetric profiles. 
We interpret v as the average rate of vertical growth due to photic response of 
microbes and A as the average rate of surface-normal growth due to mineral 
accretion. 



3 Results and discussion 



Although non- linear, equation (1) can be solved with a change of variables 
using the method of characteristics [21] and prescribed initial profiles. The 
choice of initial profile is important. Cone-like initial profiles arise naturally in 
deformations of thin fiat sheets [22]. Fig. 1 shows examples of forms obtained 
from our model using initial profiles similar to those thought in field and lab- 
oratory studies to initiate coniform stromatolites [23,24,9,25]. The functional 




(1) 
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form of these solutions to equation (1) is given by 



h{x, t) 



< 



hi{x,t), 
l + {v + Xt) 



2Xt 




Xt \p- 



I - 7 < X < -Xt\P+ 
At 1/3+1 <x< Xt\f3+ 



(2) 



hi{-x,t), 



Xt 1/3+1 < a; < At +7, 



where 



hi{x,t) 



l + {v + Xt) + 



ax^ + |At|/3+P + |/3+ 
1 - 2aAt 



X 



(3) 



with /3+ = (a ± a^)/a, 7 = min(Q;/|a|, l/a) and > |a|. The additional 
parameters, a and a, are used to tune the concavity or convexity of the flanks 
and the sharpness of the peak of the initial shape. 

The results provide possible explanations for variations in coniform stroma- 
tolite morphogenesis. When v is smaller than or comparable to A the result 
is a domical form (Figs, la, lb and Ic). Coniform structures with thickened 
apical zones form when f >> A (Figs. Id, le and If). More angular coniform 
structures form when v »> X (Figs. Ig, Ih and li). The essential charac- 
terisitics of Conophyton [26], a columnar stromatolite composed of conical 
laminae with thickened crests (Figs. 2 and 3) are apparent in Figs. Id, le 
and If. The laminae of Conophyton generally lack any evidence of trapping 
or binding of detrital sediment particles and the various microstructures that 
have been described all appear to result from a combination of BMC growth 
and carbonate precipitation [27]. The form is recorded from Paleoproterozoic 
to Mezoproterozoic rocks world-wide, but becomes rare in the Neoproterozoic 
[10]. It flourished in extensive flelds of conical columns up to 10m high in 
environments characterised by very low sedimentation rates [28,24,9,25]. The 
lack of evidence for signiflcant sedimentation and evidence for an almost com- 
plete covering by a BMC during growth suggests that the Conophyton form 
is determined by two factors, light and mineral accretion. The thickening of 
the crestal zones in Figs. Id, le and If is evocative of the thickened laminae 
and fenestrae in the delicate crestal zone of Conophyton [9]. Stromatolite fen- 
estrae are voids in the lithifled structure thought to have been left after the 
decay of the original BMC [29] . A modern analogue for Conophyton has been 
recognised in hot-springs in Yellowstone National Park, USA [23] , where it has 
been established that they form as a result of upward growth and motility of 
phototactic fllamentous cyanobacteria combined with precipitation of silica, 
and that crestal fenestrae and thickening of laminae have been related to the 
preferential upward growth of the constructing microbes. 
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It has been concluded [30] that crcstal thickening of laminae and amplification 
of bedding irregularities are evidence for phototropic growth in stromatolites. 
The results of our model, equation (1), shown in Fig. 1, together with field 
evidence, support the interpretation that the vertical growth parameter v rep- 
resents photic response of the BMC rather than sediment deposition. If the 
converse were true, coniform stromatolites would only form under conditions 
of high sedimentation which is precisely contrary to field evidence. Indeed, 
while sediment deposition would tend towards the smoothing of surface ir- 
regularities, growth due to photic response would tend to accentuate them. 
Our model shows that a combination of vertical phototropic or phototactic 
microbial growth and surface-normal mineral accretion can produce coniform 
forms and structures analogous to those found in both Archaean and Pro- 
terozoic coniform stromatolites. For example, there is a striking similarity 
between the model forms shown in Figs. Ig, Ih and li and the sharply-peaked 
coniform stromatolites in the Warrawoona Group [5], thus supporting their 
biogenic origin and reinforcing the probability that photosynthetic microbes 
were components of Archaean BMCs. The various cases modelled in Figs. Id, 
le and If can all be matched in Proterozoic Conophytons (Figs. 2 and 3). 
This sheds some light on why, after flourishing for much of the Proterozoic, 
Conophytons virtually disappeared in the Neoproterozoic [27]. This demise 
has been linked to evolutionary changes in BMCs [10], but since these would 
not have limited photic response, this seems untenable. Conophytons repre- 
sent an effective growth strategy that is especially vulnerable to predation 
and competition [14] and their demise is best explained as an indication of 
the evolution of greater biological diversity in the quiet marine environments 
that they had dominated for so long. 
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Fig. 1. Equal time snapshots of the surface height profile from the solution of the 
model, equations (1) and (2), for different values of the initial shape parameters a 
and a and the growth parameters A and First column: a = 1.5, a = 2, second 
column: a = 0.01,0; = 2, third column: a = —1.5, a = 2. First row: A = 0.5, v = 1, 
second row: A = 0.08, v = 1.5, third row: A = 0.01, v = 1.5. 
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Fig. 2. Conophytons from the 1.7 Ga Dungaminnie Formation, NT, Australia (grid 
ref. 53K NB 0578973 8154545). Coin 28mm for scale. Compare the variation in 
conical laminations with Fig. Id, le and If. 
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Fig. 3. Conophyton from the same locality as Fig. 2. Coin 28mm for scale. Note 
fenestrae, thought to have originally contained unlithified BMCs, and thickening of 
laminations in axial zone. 
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